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摘 要 核 主 泵 内 的 流动 不 稳定 将 会 引起 严重 振动 ,， 不 利于 其 安全 稳定 运行 。 因此 本 文 基于 大 涡 模 拟 (LES) 数值 计算 方 
法 对 几 个 典型 工 况 下 核 主 泵 内 部 非 稳 态 流动 结构 及 其 压力 脉动 特性 进行 全面 前述 与 关联 性 分 析 . 研究 表明 , 随 着 流量 的 增 
加 , 动静 干涉 作用 在 导 叶 出 口 处 逐渐 增加 ; 偏 工 况 条 件 下 ， 导 叶 出 口 处 压力 脉动 频谱 低频 段 中 均 出 现 复 杂 激 励 频率 ， 尤 其 
是 靠近 出 液 管 附近 的 导 叶 出 口 处 , 核 主 泵 在 但 大 流量 工 况 下 运行 时 壳 体 右 侧 内 部 非 稳 态 流动 结构 相 较 于 这 体 左 侧 更 加 复 
杂 ; 在 仿 小 流量 工 况 下 运行 时 壳 体 底部 压力 脉动 更 加 剧烈 , 本 文 进 一 步 详细 描述 了 核 主 泵 球形 过 体内 强 涡 量 区 的 流动 结构 
及 其 成 因 , 并 且 发 现 测 点 处 的 压力 频谱 与 涡 量 频谱 有 相同 的 主要 激励 频率 ,因此 证 明 核 主 泵 内 非 定常 旋涡 流动 结构 是 激励 
压力 脉动 的 主因 之 一 。 
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Abstract Flow instability will cause severe vibration in the nuclear reactor coolant pump (RCP), 
and it is bad for the safe and stable operation. Thus, this article is based on large eddy simulation 
(LES) numerical calculation method to study the internal unsteady flow structures and pressure 
pulsation characteristics and correlation analysis of the RCP under several typical working conditions. 
Research shows that with the increase of flow rate, the rotor-stator interactions gradually increase 
in the diffuser outlet. Under the partial working conditions, several complex frequencies are found 
in low frequency band of pressure spectrum, especially the diffuser outlet near the discharge nozzle. 
The internal unsteady flow structures in the right side of spherical casing are more complicated than 
the left side under the large flow conditions. The pressure fluctuation is more intense at the bottom 
of the spherical casing under the small flow conditions. Further, the unsteady flow structures and 
its reasons in high vorticity magnitude region are detailly described in the ROP. It has been found 
that the vorticity spectra and the pressure spectra almost have the same main excitation frequencies. 
Therefore, it has been confirmed that for particular regions pressure pulsations are determined by the 
vortex shedding of the RCP. 
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0 引言 


由 于 国外 对 核反应 堆 冷 却 剂 循环 泵 (以 下 简称 
核 主 泵 ) 的 技术 垄断 , 我 国 在 研制 高 效 核 主 泵 技术 方 
面 成 为 了 核电 技术 产业 化 的 瓶 宽 . 因此 , 研究 高 效 且 
稳定 运行 的 核 主 泵 不 仅 具 有 重要 的 节能 意义 ,同时 
对 我 国 目 主 研发 核电 技术 也 具有 更 加 深 二 的 影响 ， 
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从 基础 性 研究 角度 出 发 ， 铬 要 设计 高 效 稳定 的 核 主 
泵 ， 就 必须 从 最 基础 的 泵 内 部 非 稳 态 流 动 入 手 ， 还 
原核 主 泵 内 部 复杂 流动 结构 ,掌握 非 稳 态 流动 结构 
与 激励 特性 之 间 的 关联 关系 ， 同 时 寻找 流动 损失 的 
主要 来 源 也 是 核 主 泵 研究 的 必然 趋势 。 
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核 主 泵 内 部 流动 不 稳定 特性 会 导致 严重 的 振 
动 是 ,尤其 在 特殊 的 准 球 型 压 水 室 中 外 ,内 部 非 稳 
态 流 场 及 其 复杂 ， 并 且 对 其 稳定 运行 产生 危害 。 相 
较 于 普通 螺旋 形 压 水 室 8 , 准 球形 压 水 室内 会 产生 
典型 的 流动 分 离 、 二 次 流 及 回流 现象 , 并 且 压 水 室 两 
侧 的 非 定 常 流动 结构 极其 不 对 称 可 .Claus Knierim 
等 口 为 1400 MW 的 核电 站 (APR1400) 自主 研发 了 
一 台新 型 核 主祭, 利用 计算 流体 动力 学 (CFD) 技术 
优化 了 叶轮 与 导 叶 ， 使 其 性 能 达到 最 优 .在 研究 过 
程 中 ,他 们 发 现 核 主 泵 内 的 流体 堵塞 在 位 于 出 液 管 
附近 的 导 叶 流 道 中 有 旦 流体 流动 发 生 分 离 现象 。 其 中 
一 部 分 流体 流出 出 液 管 ， 另 一 部 分 流体 在 球形 压 水 
室内 再 次 循环 ， 这 也 是 造成 壳 体 内 损失 的 主要 因素 
之 一 回 . Sven Baumgarten 等 上 详细 介绍 了 一 台 用 
于 西屋 AP1000 反应 堆 的 RUV 型 核 主 泵 的 结构 . 其 
水 力 部 件 的 设计 是 基于 CFD 技术 完成 的 , 并 且 最 终 
采用 1:2 的 模型 泵 展开 了 全 面 的 试验 研究 ,韩国 原 
子 能 研究 院 (KAERI) 为 1400 MW 核反应 堆 设计 
制造 了 一 台 测 试用 核 主 泵 ,进行 了 冷 态 与 热 态 性 能 
测试 ， 观 察 到 在 特定 的 温度 范围 内 ， 压 力 脉动 振幅 
异常 增加 。 其 原因 是 当 核 主 泵 叶片 通过 频率 及 其 谐 
波 与 核 主 泵 测试 设备 的 频率 成 正比 时 ， 表 共振 现象 
发 生 . Kato 等 8 站 采用 大 涡 模 拟 方 法 描述 了 在 小 流 
量 工 况 下 的 混流 泵 的 内 部 流 场 。 结 果 表 明 预 测 流 量 
-扬程 曲线 与 测量 结果 相符 合 ， 其 扬程 - 流量 曲线 
在 小 流量 工 况 下 呈现 不 稳定 特性 . Miyabe 等 0 
运用 PIV 和 压力 脉动 测量 仪器 对 混流 泵 内 部 旋转 失 
速 现象 进行 了 研究。 研究 表明 不 稳定 内 部 流动 特性 
是 由 在 叶轮 出 口 至 导 叶 进口 处 出 现 的 大 尺度 回流 所 
激励 的 。 以 上 研究 显示 ， 核 主 泵 内 非 定 党 流动 结构 
可 能 会 导致 严重 的 压力 脉动 和 振动 ， 将 会 危害 泵 的 
安全 稳定 运行 吧 9 

由 于 核 主 泵 球形 党 体 结构 的 特殊 性 ， 其 内 部 的 
不 稳定 流动 结构 与 激励 特性 之 间 的 关联 关系 尚未 被 
完全 揭示 ， 因 此 本 文 以 已 有 的 核 主 泵 模型 为 研究 对 
和 象 , 运用 大 涡 模拟 (LES) 方法 模拟 了 在 几 个 典型 运 
行 工 况 下 核 主 泵 内 非 定常 流动 特性 与 其 激励 特性 ， 
为 核 主 泵 安全 稳定 运行 英 定 基础 . 


1 模型 到 主要 参数 及 计算 区 域 


1.1 模型 条 主要 参数 

以 已 有 的 混流 式 核 主 泵 模型 为 研究 对 每 ， 其 主 
要 设计 参数 和 几何 参数 如 表 1 所 示 . 核 主 泵 的 水 力 
部 件 主要 有 叶轮 、 导 叶 和 球 元 组 成 。 为 满足 主 泵 中 
老 体 承 压 安全 性 能 好 、 强 度 高 、 易 加 工 等 要 求 , 采用 


677 
准 球形 壳 休 的 设计 . 
表 1 模型 泵 主要 参数 


Table 1 Designing parameters 


参数 表达 式 与 数值 
额定 流量 Qa/(m3/h) 848 
设计 扬程 Ha/m 12.7 
额定 流量 系数 Ba Qa/u2 R2=0.63 
额定 扬程 系数 wa gHa/u3=0.29 
额定 转速 ng/(r/min) 1480 
叶轮 叶片 数 孔 4 
导 叶 叶片 数 2a 12 
叶轮 进口 直径 D1/mm 221 
叶轮 出 口 直径 D2/mm 268 
叶轮 出 口 宽 度 bz /mm 84 
球形 壳 体 直径 D3/mm 637.5 


1.2 计算 区 域 

图 1 给 出 了 模型 泵 的 计算 区 域 , 静止 域 分 别 为 
进口 流域 、 导 叶 流 道 、 球形 党 体 流 域 和 出 口 流域 , 转 
动 域 为 叶轮 流 道 . 


出 口 


| 球形 壳 体 


图 1 计算 区 域 示意 图 


Fig. 1 Three-dimensional computational domain 


2 数值 计算 
2.1 数值 计算 方法 

本 文采 用 LES 数值 计算 方法 对 核 主 泵 在 不 同 
工 况 下 的 非 定常 流动 结构 及 其 激励 特性 进行 计算 分 
析 . LES 方法 的 基本 思路 是 : 用 瞬时 Navier-Stokes 
方程 直接 模拟 灌流 中 的 大 尺度 涡 ， 不 直接 模拟 小 尺 
度 涡 , 然后 通过 求解 附加 方程 得 到 小 尺度 涡 的 解 . 显 
然 这 种 求解 方法 相 比 于 直接 数值 模拟 效率 更 高 ， 消 
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耗资 源 更 少 , 但 却 比 RANS 方法 更 加 精确 04. 因此 
LES 数值 计算 方法 相 较 于 RANS 计算 方法 有 更 好 的 
能 力 预测 瞬 态 流 场 , 尤其 在 流动 分 离 区 域 0519. 本 
文 应 用 商业 软件 Ansys-Fluent 13.0 预测 核 主 奈 中 三 
维 非 定常 流动 特性 . 为 得 到 更 精确 计算 结果 , 首先 采 
用 标准 h-e 计算 模型 对 核 主 泵 进行 定常 计算 , 将 定 
常 计算 结果 作为 非 定常 计算 结果 的 初始 条 件 。 对 于 
非 定常 计算 ， 采 用 Smagorinsky-Lilly 亚 格 子 (SGS) 
混合 系数 模型 047 对 方程 进行 封闭 . 
2.2 网 格 划 分 与 边界 条 件 

数值 计算 的 精度 主要 取决 于 网 格 质 量 . 全 局 网 
格 的 生成 应 用 结构 化 网 格 技术 (ICEM 网 格 生 成 工 
具 ) 进行 网 格 划 分 ， 并 且 在 叶片 表面 等 局 部 位 置 进 
行 加 密 处 理 。 因 此 将 获得 更 加 精确 的 核 主 泵 内 部 
非 稳 态 流动 结构 与 压力 脉动 特性 .同时 为 兼顾 计算 
准确 性 以 及 计算 机 的 计算 能 力 ， 通 过 网 格 无 关 性 验 
证 ， 最 终 确 定 全 局 网 格 节 点 数 为 7163843， 网 格 数 
为 7438559， 其 中 时 轮 节 点 数 为 3127239， 网 格 数 为 
3231400. 图 2(a) 为 整 泵 网 格 图 , 图 2(b) 为 叶轮 网 格 
及 局 部 加 蜜 网 格 示 意图 . 

进出 口 分 别 定义 速度 进口 边界 条 件 (velocity- 
inlet) 和 压力 出 口 边界 条 件 (pressure-outlet)， 压 力 
为 环境 压力 . 非 定 常 计 算 时 ， 旋 转 域 (叶轮 ) 与 静止 
域 ( 导 叶 ) 之 间 的 耦合 采用 滑 移 网 格 模型 ， 进 出 口 流 
道 、 导 叶 流 道 和 球 壳 流域 均 和 采用 静止 壁面 无 滑 移 边 
界 条 件 , 在 近 壁 区 采用 标准 壁面 孔 数 法 ， 


(b) 叶轮 网 格 及 局 部 加 密 


图 2 ”网 格 示意 图 
Fig. 2 Grids 


2.3 非 定常 设置 

为 精确 辨析 出 球形 元 体内 流动 的 非 定 稼 特性 ， 
本 文 将 时 间 步 长 设置 为 At=112.6 us， 即 每 一 时 间 
步 长 叶轮 转 过 的 角度 为 *，360 个 步 长 叶轮 转动 一 
周 . 同时 计算 残 差 收敛 性 依据 设 定 为 9x10 来 获 
得 更 加 准确 的 计算 结果 。 另外， 设置 叶轮 旋转 五 十 
圈 后 可 认定 为 稳定 可 靠 的 计算 结果 . 


图 3 监测 点 分 布 


Fig. 3 Monitoring points on each main flow domain 


2.4 监测 点 分 布 

为 了 全 面 获得 球形 元 体内 的 压力 脉动 特性 ,在 
核 主 泵 中 间 规 面 的 重要 位 置 设置 压力 脉动 监测 点 ， 
图 3 给 出 了 监测 点 位 置 布置 图 . 球形 元 体 的 主要 作 
用 为 承 压 , 其 周 癌 的 压力 脉动 分 布 极为 重要 , 因此 在 
元 体 周 癌 每 45° 设 一 个 监测 点 , 分 别 为 监测 点 W1- 
W7. 导 叶 尾 缘 容 易 出 现 周 期 性 旋涡 脱落 ,其 脱落 特 
性 与 压力 脉动 特性 关联 极 大 ， 因 此 在 导 叶 尾 缘 处 每 
90° 设 一 个 监测 点 , 分 别 为 监测 点 D1-D4. 


3 计算 结果 分 析 
3.1 非 定常 压力 脉动 分 析 


为 了 分 析 非 定常 压力 脉动 变化 情况 ， 引 入 压力 
系数 Cp, 党 用 的 压力 系数 表达 式 为 : 


AP 


Cp= 卫 (1) 
二 Du22 
2 
7T Do 
2 2 
U2 60 (2) 


式 中 ，AP 为 监测 静 压 与 其 平均 压力 之 差 ; va 为 叶 
轮 出 口 圆周 速度 ; Ds 为 叶轮 出 口 直径， 

叶轮 转速 为 1480 r/min,， 所 以 轴 频 ( 户 ) 与 叶轮 
通过 频率 (fepF) 分 别 为 24.7 Hz 和 99 Hz。 本文 运 
用 快速 傅 里 时 (FFT) 方法 将 时 域 信号 转换 为 频 域 信 
号 。 在 信号 处 理 过 程 中 ,采用 汉 宁 徐 (Hanning) 来 
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图 4 各 工 况 下 导 叶 出 口 处 压力 脉动 频谱 


Fig. 4 Pressure pulsation spectra on diffuser outlet monitoring points 


减少 能 量 泄漏 .为 了 全 面 了 解 导 叶 出 口 压力 脉动 特 
性 , 图 4 给 出 了 导 叶 出 口 处 压力 脉动 频谱 分 布 。 从 
图 4 可 以 看 出 ,在 额定 工 况 (8Q/Qa=1.0) 下 ， 各 测 
点 的 主要 激励 频率 均 为 一 倍 叶 频 (fppr) 与 二 倍 叶 频 
(2 户 PF). 因此 在 额定 工 况 下 导 叶 出 口 处 影响 压力 脉 
动 的 主要 因素 是 动静 干涉 作用 .值得 注意 的 是 ,在 
D1 测 点 靠近 出 液 管 处 的 导 叶 出 口 处 在 低频 段 出 现 
一 些 激励 频率 ， 这 说 明 在 额定 工 况 下 D1 测 点 导 叶 
出 口 处 压力 脉动 相 较 于 其 他 测 点 更 加 复杂 . 在 小 流 
量 工 况 下 (8Q/Qa=0.7)， 激 励 频率 除了 叶 频 (JPF) 
外 ， 在 低频 段 均 出 现 了 主要 激励 频率 ， 其 频率 区 间 
约 在 8~50 Hz, 这 说 明 偏 小 流量 工 况 时 , 导 叶 出 口 处 
的 非 定 常 流动 结构 所 引起 的 压力 脉动 抑制 了 动静 干 
涉 作 用 . 在 大 流量 工 况 下 (Q/Qa=1.2), 各 测 点 在 叶 
频 (fepF) 处 ,压力 脉动 幅 值 较 其 他 工 况 均 增加 , 由 
此 说 明 偏 大 流量 工 况 下 导 叶 出 口 处 动静 干涉 作用 更 
加 明显 . 并且 在 低频 段 8~30 Hz 频率 区 间 内 ,出 现 
了 较 多 高 幅 值 的 非 定 常 激 励 频 率 , 尤其 在 监测 点 D1 
处 非 定 常 激励 频率 更 加 复杂 。 综 上 所 述 ， 随 着 流量 
的 增加 ， 动 静 干 涉 作 用 在 导 叶 出 口 处 逐渐 增加 。 偏 
工 况 条 件 下 ， 低 频段 均 出 现 复杂 激励 频率 ， 尤 其 是 
靠近 出 液 管 附近 的 导 叶 出 口 处 ， 压 力 脉动 相 较 其 他 
测 点 处 更 加 复杂 ， 

为 深入 分 析 壳 体 周 向 压力 脉动 分 布 特性 , 图 5 
与 图 6 分 别 给 出 了 各 工 况 下 壳 体 周 向 各 测 点 一 倍 
叶 频 处 压力 脉动 幅 值 分 布 与 压力 脉动 幅 值 均 方 根 误 
差 (RMSE) 分 布 。 其 中 压力 脉动 幅 值 用 均 方 根 误差 
(RMSE) 来 评价 , 其 公式 为 


-0/04=0.7 
一 -0O10i=1.0 
— 和 —O/O=l1 2 


0.004 Wl 


图 5 壳 体 周 向 各 测 点 一 倍 叶 频 处 压力 脉动 幅 值 分 布 
Fig. 5 Pressure pulsation amplitudes offBpF at monitoring 
points Wl-W7 


其 中 ,Cb 是 压力 系数 平均 值 , ”是 所 测量 压力 系数 
的 数量 . 因此 均 方 根 误差 (RMSE) 代表 压力 脉动 相 
较 于 平均 值 的 离散 程度 。 从 图 5 可 以 看 出 ,各 工 况 
下 监测 点 在 圆周 方向 上 的 动静 干涉 作用 (fppr) 较为 
不 对 称 。 动静 干涉 作用 (fepr) 在 测 点 W4 处 幅 人 
均 达 到 最 大 ， 尤 其 是 在 大 流量 工 况 下 (Q/Qa=1.2)， 
这 说 明 核 主 泵 在 偏 大 流量 工 况 下 运行 时 动静 干涉 对 
球形 元 体 底部 作用 更 加 明显 . 此 外 ， 从 图 6 元 体 周 
加 各 测 点 压力 脉动 幅 值 均 方 根 误差 分 布 可 以 看 出 ， 
额定 工 况 与 俩 大 流量 工 况 下 压力 脉动 幅 值 在 元 体 周 
加 分布 情况 相 一 致 ， 但 小 流量 工 况 下 脉动 幅 值 分 布 
差别 较 大 . 人 额定 工 帝 与 偏 大 流量 工 况 下 压力 脉动 幅 
值 在 壳 体 右 侧 测 点 W6 处 达到 最 大 ， 党 体 右 侧 压 
力 脉动 幅 值 均 大 于 党 体 左 侧 ， 这 说 明 核 主 泵 在 偏 大 


201710.00081v1 


chinaXiv 


1680 工程 热 物 理 学 报 


流量 工 况 下 运行 时 壳 体 右 侧 内 部 非 稳 态 流动 相 较 于 


中 间 截 面 的 瞬 态 涡 量 分 布 .将 强 涡 量 区 分 为 a 与 


元 体 左 侧 更 加 复杂 。 在 小 流量 工 涡 下 ， 各 测 点 压 
力 脉动 幅 值 除 W4 测 点 外 基本 对 称 分 布 。W4 测 
点 脉动 幅 值 在 小 流量 工 况 下 达到 最 大 ， 这 说 明 核 主 
泵 在 偏 小 流量 工 况 下 运行 时 腕 体 底部 压力 脉动 更 加 
剧烈 ， 


Cn 
0.012 
一 0/0r07 


-全 _O/1OD=1.2 
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图 6 壳 体 周 向 各 测 点 压力 脉动 幅 值 均 方 根 误差 分 布 
Fig. 6 RMSE values of monitoring point VV1-VV7 


3.2 内 流 结构 与 压力 脉动 特性 关联 关系 
为 了 分 析 核 主 奈 内 部 非 稳 态 流 动 结构 与 压力 脉 
动 之 间 的 关联 关系 , 图 7 给 出 了 三 个 工 况 下 核 主 泵 


两 个 区 域 ， 并 给 出 了 这 两 个 区 域 的 局 部 放大 图 .出 
前 文 的 分 析 可 知 ， 导 叶 叶 片 A 尾 缘 D1 测 点 处 压力 
脉动 在 低频 段 较为 复杂 ,同时 学 体 右 侧 压力 脉动 幅 
值 较 左 侧 更 大 ， 因 此 将 接近 出 液 管 附近 的 导 叶 叶片 
分 为 a 区 ; 在 偏 大 流量 工 况 下 , 党 体 右 侧 W6 测 点 
处 压力 脉动 幅 值 最 大 ， 因 此 将 右 侧 强 涡 量 区 划分 为 
6 区 域 . 

在 额定 工 况 下 (1.084), 从 a 区 可 以 看 出 , 在 导 
叶 叶 片 A 与 B 吸力 面 中 缘 至 尾 缘 处 出 现 了 明显 的 
非 稳 态 旋涡 区 ， 并且 在 尾 缘 出 现 了 旋涡 脱落 。 这 是 
因为 导 叶 叶片 A 与 B 正 对 着 出 液 管 , 其 液 流 流出 导 
叶 后 由 于 惯性 的 作用 撞击 到 球形 腕 体 右 侧 出 液 管 壁 
面 ,形成 了 大 太 度 的 流动 分 离 ， 其 中 一 部 分 液 流 随 
着 主流 方向 ( 顺 时 针 方 向 ) 流 去 , 而 男 一 部 分 流体 在 
压力 梯度 与 堪 侧 流体 的 挤 压 作用 下 流出 出 液 管 . 在 
此 撞击 与 挤 压 作用 下 ， 致 使 导 叶 叶片 A 与 B 尾 缘 
形成 非 稳 态 旋涡 结构 ， 因 此 在 D1 测 点 处 产生 了 较 
大 的 压力 脉动 . 在 9 区 中 , 在 叶轮 顺 时 针 的 扫 掠 下 ， 
叶轮 出 口 处 形成 了 明显 的 射流 尾 迹 ， 其 油 量 相对 较 
大 . 在 叶轮 射流 尾 迹 与 导 叶 的 动静 干涉 下 ， 其 尾 迹 
不 断 通 过 导 叶 各 流 道 , 在 导 叶 各 流 道中 形成 了 涡 市 ， 
此 涡 带 随 着 主流 移动 到 导 叶 出 口 逐渐 扩散 ， 并 在 导 
叶 叶 片尾 缘 形 成 了 类 似 卡门 涡 街 的 周期 性 旋涡 脱落 
结构 ， 


图 7 各 工 沉 下 核 主 泵 中 间 截 面 的 瞬 态 涡 基 分 布 


Fig. 7 Transient vorticity distributions of three conditions in the middle section 
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图 8 上 额定 工 况 下 D1 测 点 的 压力 频谱 ( 左 图 ) 与 涡 量 频谱 ( 右 图 ) 
Fig. 8 Pressure spectrum (left) and vorticity spectrum (right) of point D1 under the rated conditions 
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图 9 额定 工 况 下 D2 测 点 的 压力 频谱 ( 左 图 ) 与 涡 量 频谱 ( 右 图 ) 


Fig. 9 Pressure spectrum (left) and vorticity Spectrum (right) of point D2 under the rated conditions 


在 大 流量 工 况 下 (1.2Qa), 叶轮 射流 尾 迹 所 形成 
的 涡 量 相 较 额定 工 况 进一步 增加 ， 射 流 尾 迹 结 构 在 
导 叶 流 道内 形成 的 涡 带 变 宽 。 从 a 区 中 可 以 看 出 ， 
导 叶 叶片 A 和 B 的 吸力 面 中 部 至 尾 缘 处 的 强 滴 量 
区 进一步 扩大 , 非 稳 态 旋 涡 结构 扩散 至 出 液 管 右 侧 ， 
这 些 非 稳 态 旋 涡 结构 将 会 堵塞 导 叶 叶片 A 与 B 之 间 
的 流 道 , 甚至 阻碍 流体 流出 出 液 管 . 除了 接近 出 液 管 
附近 的 导 叶 叶片 A 和 B 中 部 至 尾 缘 处 的 非 定常 旋涡 
结构 ， 其 他 导 叶 叶片 吸力 面 中 部 至 尾 缘 处 均 出 现 强 
涡 量 区 . 从 0 区 中 可 以 看 出 , 大 流量 工 况 下 导 叶 叶 
片 D 吸力 面 中 部 附近 边界 层 开 始 分 离 , 并 且 形 成 旋 
涡 区 , 旋涡 区 一 直 延 伸 至 导 叶 叶片 尾 缘 ， 直 至 脱落 . 
大 流量 工 况 下 ， 尾 缘 周 期 性 旋涡 脱落 结构 相 较 于 额 
定 工 况 更 加 明显 ， 

在 小 流量 工 况 下 (0.784a), 核 主 泵 内 部 非 稳 态 流 
动 结构 相 较 于 其 他 工 况 更 加 复杂 . 总 体 上 ， 壳 体内 
涡 量 强度 远大 于 其 他 两 工 况 ， 并 且 所 有 导 叶 流 道 均 
堵塞 严重 。 在 叶轮 叶片 出 口 边 靠近 工作 面 处 ， 相 对 
流速 大 于 靠近 叶片 背面 的 流体 速度 ， 因 此 在 叶片 出 
口 处 形成 了 较 其 他 两 个 工 况 更 加 明显 的 非 均 匀 射 流 
尾 迹 结构 。 这 样 强 涡 量 的 射流 尾 迹 结构 随 着 动静 干 
涉 作用 下 流入 各 导 叶 流 道 ， 使 得 导 叶 各 流 道中 的 流 
体 相互 挤 压 堵 塞 ,使 得 流体 不 能 顺利 通过 导 叶 流 道 ， 
这 些 复杂 非 定常 涡 旋 结构 在 惯性 力作 用 下 逐渐 向 导 
叶 出 口 处 扩散 , 引起 了 导 叶 出 口 处 强烈 的 压力 脉动 ， 
造成 不 同 测 点 处 压力 信号 在 低频 段 8~50 Hz 处 出 现 


了 较 多 高 幅 值 的 复杂 激励 频率 ， 

为 了 进一步 分 析 非 稳 态 流动 结构 与 压力 脉动 之 
间 的 关系 ， 以 额定 工 况 为 例 , 图 8 给 出 了 该 工 况 下 
D1 测 点 的 压力 频谱 ( 左 图 ) 与 涡 量 频谱 ( 右 图 ). 从 图 
中 可 以 发 现 , 在 D1 测 点 压力 频谱 与 涡 量 频谱 有 相同 
的 主要 激励 频率 , 除了 叶 频 外 其 主要 激励 均 为 25 Hz 
与 32 Hz. 由 此 可 以 进一步 证 明 泵 内 非 定常 周期 性 旋 
涡流 动 结构 是 激励 压力 脉动 的 主要 因素 之 一 . 图 9 
给 出 了 额定 工 况 下 D2 测 点 的 压力 频谱 ( 左 图 ) 与 涡 
量 频 谱 ( 右 图 ). 由 于 D2 测 点 位 于 导 叶 叶片 D 出 口 
周期 性 旋涡 脱落 处 , 此 时 D2 测 点 压力 频谱 与 涡 量 频 
谱 的 主要 激励 频率 均 为 1 倍 叶 频 (fppr), 尤其 在 涡 
量 频 谱 中 , 1 倍 叶 频 处 脉动 幅 值 较 大 . 由 此 现象 可 认 
为 ， 导 叶 尾 缘 周 期 性 旋涡 脱落 频率 为 1 倍 叶 频 . 这 
是 因为 导 叶 流 道中 的 涡 带 在 动静 干涉 作用 下 疝 导 叶 
出 口 不 断 推进 ， 当 推进 至 导 叶 叶片 尾 缘 人 处 与 尾 缘 处 
旋涡 相合 并 , 随 之 脱落 . 


4 结 论 


本 文采 用 LES 数值 计算 方法 对 核 主 泵 在 不 同 工 
况 下 的 非 定 常 流动 结构 及 其 激励 特性 进行 计算 分 析 ， 
得 到 以 下 结论 : 

1) 随 着 流量 的 增加 ,动静 干涉 作用 在 导 叶 出 口 
处 逐渐 增加 。 俩 工 况 条 件 下 ， 低 频段 均 出 现 复杂 激 
励 频 率 ， 尤其 是 靠近 出 液 管 附近 的 导 叶 出 口 处 ， 压 
力 脉动 相 较 其 他 测 点 处 更 加 复杂 ; 
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2) 核 主 泵 在 偏 大 流量 工 况 下 运行 时 壳 体 右 侧 内 
部 非 稳 态 流动 相 较 于 元 体 左 侧 更 加 复杂 ; 在 侦 小 流 
量 工 况 下 运行 时 腕 体 底部 压力 脉动 更 加 剧烈 ; 

3) 在 额定 工 况 下 (1.08a), 在 导 叶 叶片 A 与 B 
吸力 面 中 部 至 尾 缘 处 出 现 了 明显 的 非 稳 态 旋 涡 结构 ， 
并 且 在 尾 缘 出 现 了 旋涡 脱落 。 在 叶轮 射流 尾 迹 与 导 
叶 的 动静 干涉 下 ， 其 尾 迹 不 断 通 过 导 叶 各 流 道 ,在 
导 叶 各 流 道中 形成 了 涡 带 ,这 样 的 涡 带 随 着 主流 移 
动 到 导 叶 出 口 逐渐 扩散 ,并 在 导 叶 叶片 尾 缘 形 成 了 
周期 性 旋涡 脱落 结构 ; 

4) 在 大 流量 工 沉 下 (1.28Qa), 叶轮 射流 尾 迹 所 形 
成 的 涡 量 相 较 额定 工 况 进一步 增加 ， 射流 尾 迹 结构 
在 导 叶 流 道内 形成 的 镜 囊 变 宽 ， 导 叶 尾 绿 周期 性 旋 
涡 脱落 结构 相 较 于 额定 工 帝 更 加 明显 ; 

5) 在 小 流量 工 况 下 (0.784)， 核 主 泵 内 部 非 稳 
态 流动 结构 相 较 于 其 他 工 况 更 加 复杂 ， 沉 体内 涡 量 
强度 远大 于 其 他 两 工 况 ， 并 且 所 有 导 叶 流 道 均 堵 赛 
严重 ; 

6) 测 点 处 的 压力 频谱 与 涡 量 频谱 有 相同 的 主要 
激励 频率 ,因此 泵 内 非 定常 周期 性 旋涡 流动 结构 是 
激励 压力 脉动 的 主要 因素 之 一 ， 
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